The mystery of O and O3 production in the effluent of a 
He/02 atmospheric pressure microplasma jet 

D EUerweg, A von Keudell, and J Benedikt 

Institute for Experimental Physics II: Reactive Plasmas, 
Rxihr-Universitat Bochum, 44780 Bochimi, Germany 

E-mail: jan.benedikt@ruhr-uni-bochuin.de 

Abstract. Microplasma jets are commonly used to treat samples in ambient air 
atmosphere. The effect of admixing air into the effluent may severely affect the 
composition of the emerging species. Here, the effluent of a He/02 microplasma jet 
has been analyzed in a helium and in an air atmosphere by molecular beam mass 
spectrometry. First, the composition of the effluent in air has been recorded as a 
frmction of the distance to determine how fast air admixes into the effluent. Then, the 
spatial distribution of atomic oxygen and ozone in the effluent has been recorded in 
ambient air and compared to measurements in a helium atmosphere. Additionally, 
a fluid model of the gas flow with reaction kinetics of reactive oxygen species in the 
effluent has been constructed. In ambient air, the O density declines only slightly 
faster with the distance compared to a helium atmosphere. On the contrary, the O3 
density in ambient air increases significantly faster with the distance compared to a 
helirmi atmosphere. This mysterious behavior can have big implication for the use of 
similar jets in plasma medicine. It is shown that photodissociation of 02 and 03 is 
not responsible for the observed effect. A reaction scheme involving the reaction of 
plasma produced highly vibrationally excited O2 with ground state O2 molecules is 
proposed as a possible explanation of the observed densities. A very good agreement 
between measured and simulated densities is achieved. 
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1. Introduction 

Cold atmospheric plasmas exhibit many unique properties, which make them very 
attractive for a broad field of applications. These plasmas offer non-equilibrium 
chemistry at atmospheric pressure with high densities of reactive neutral and charged 
species and high fluxes of energetic photons. They can be utilized as light sources, 
for localized treatment of thermo and vacuum sensitive materials, and, especially in 
the case of jets, also of living tissue in plasma medicine applications. 
Cold atmospheric plasma jets have been developed and studied in the past for many 
different applications, like thin film deposition (H. 12113]], inactivation of pathogens 
im ]5l ]6l ]Z]] or generation of nanoparticles |]8]]. Many of the developed sources use 
helium as plasma forming gas |]9l]10l]lll]12l with a small admixture of a molecular gas 
(e.g. molecular oxygen). It has been shown that these jets are effective in producing 
reactive oxygen species (ROS) and can effectively inactivate for example bacteria or 
etch hydrocarbon polymer layers 113 . 

For the understanding of the interaction of the plasma effluent with substrates, of 
inactivation of bacteria, or of its effects on living tissues, it is necessary to know 
absolute densities and fluxes of reactive species in the plasma jet. Moreover, the 
microplasma jets are in many cases operated in ambient atmosphere, which effect on 
the plasma effluent is not fully and quantitatively understood. 

Measurements of absolute densities of reactive species in microplasma jets are a 
nontrivial task due to the high pressure and the small dimensions. Advanced 
diagnostics like two-photon absorption laser-induced fluorescence spectroscopy 
(TALIF) [)TOl ]T4l ]T5| , molecular beam mass spectrometry (MBMS) p5| or absorption 
spectroscopy IIT6I have to be used to obtain absolute ROS densities. It is important 
for possible applications that all measurements are carried in realistic situations. 
Many experiments are, however, performed in a controlled atmosphere of noble gas, 
whereas possible applications will be performed in ambient air. 
Here, absolute atomic oxygen and ozone densities in the effluent of a He/02 micro- 
scaled atmospheric pressure plasma jet (/.i-APPJ) operated in ambient helium and air 
are measured by absolutely calibrated MBMS. The main focus of the work presented 
here is to analyze the influence of ambient air on the ROS densities and to reveal 
possible reaction pathways leading to their production. This is achieved by a 
comparison of the experimental data to a fluid model of the gas flow combined 
with chemical reactions in the plasma effluent. 

2, Experimental setup 

Only a brief description of the experimental setup is provided here since it has been 
described in detail elsewhere |]T5l ]T7| . The main difference is that in this study 
the jU-APPJ has been operated and analyzed in ambient air in contrast to previous 
measurements in an ambient helium atmosphere. 
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Figure 1. Sketch (a) and photo (b) of the fi-APPJ 



2.1. Micro-scaled atmospheric pressure plasma jet 

A schematic sketch of the already defined ft-APPJ is shown in figure [H The 
microplasma jet consists of two parallel stainless steel electrodes (30 mm long, 1 mm 
thick) separated by a gap of 1 mm. One electrode is powered by an rf power supply 
(13.56 MHz, absorbed power <1 W), while the other one is grounded. The /j-APPJ is 
operated at a gas flow of 1.4 slm helium with a small admixture of molecular oxygen 
(<1.6 %). The plasma volume of 30x1x1 mm^ is confined on both sides by glass plates. 
If not stated otherwise, the plasma is operated with a gas flow of 1.4 slm He with 
0.6 % O2 and an applied electrode voltage of 230 Vrms. 

Previous measurements in a He atmosphere IITOi IT4l ITSl revealed that the /.i-APPJ is 
very effective in producing ROS. An O2 dissociation rate up to 20% and maximum 
O and O3 densities around lO^^cm"^ could be measured. The O density is maximal 
at an O2 admixture of 0.6 %, while the O3 density continuously increases with the O2 
admixture. The atomic oxygen density increases when a higher power is applied to 
the jU-APPJ, whereas the ozone density slightly decreases. With increasing distance 
from the ft-APPJ nozzle, the O concentration decreases and the O3 concentration 
increases. 



The so-called X-Jet is a modified version of the /j-APPJ and allows the separation of 
heavy species from photons in the effluent downstream of the plasma. The separation 
is realized in the crossed channel structure at the nozzle of the jet, where a direct 
channel (extension of the space between electrodes) is crossed by an additional side 
channel (cf. figure IZj- A helium flow introduced through the side channel of this 
structure deflects the particles in the plasma effluent (incl. radicals and metastables) 
from their movement through the direct channel into a side channel. The photons, 
on the contrary, are not affected by this additional gas flow and they can propagate 
undisturbed through the direct channel. A detailed description of the X-Jet and tests 
of its performance can be found in the literature | [T3l[T8| . As will be demonstrated later. 



2.2. X-Jet 
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Figure 2. Sketch (a) and photo (b) of the X-Jet 



the X-Jet allows to test the photochemistry induced by plasma generated photons in 
the gas mixture introduced through the side channel (cf . figure |2]a) and in ambient 
air. 



2.3. Molecular beam mass spectrometer 

A molecular beam mass spectrometer (MBMS) is used to study the effluent of the 
/.i-APPJ. Because the operation of a mass spectrometer (MS) requires a good vacuum, 
a differential pumping system has to be used. The system has been described in detail 
previously HZl. In summary, the MBMS consists of three pumping stages that are 
connected by small orifices (cf. figure [S]). The gas extraction from the atmosphere 
into the first pumping stage is done by a sampling orifice with a diameter of 100 [,im. 
In the first stage, a special beam chopper (i.e. rotating flat metal disk with four small 
embedded skimmers) is mounted to allow only a pulsed gas flow into the differential 
pumping system. The second and third pumping stages are interconnected via a 
skimmer (0 0.8 mm). The ionizer of the MS is installed in the third stage in Itne- 
of-sight with the orifice and the skimmer to ensure that the molecular beam (MB) 
forming behind the sampling orifice can reach the ionizer without any disturbances. 
Two different MBMS designs have been used. On the one hand, the previously 



described one [TZl with the front plate of the MBMS oriented vertically like shown in 
figure |3] and on the other hand a new MBMS setup with slightly different dimension 
and a horizontally oriented front plate. The axis of the /,t-APPJ has always been placed 
normal to the sampling orifice pointing directly at it (see figure ID). 
A chamber for a controlled atmosphere is mounted at the front plate of the MBMS 
system and encloses the plasma jet. The volume of the chamber (approx. 1.1 dm'^) 
can be filled with helium or air. 

The MS measures the species density in the ionizer. With a proper background 
correction, the composition of the MB in the ionizer can be determined. When gas is 
sampled from atmosphere to low pressure, the composition of the MB is, however, 
different from the original composition of the gas mixture at the sampling orifice. 
This is due to the so-called composition distortion in MBMS sampling. A supersonic 



O and O3 production in the effluent of a He/02 microplasma jet 



5 



rotating 
skimmer 




Figure 3. Schematic of the differential pumping system 



expansion takes place downstream from the sampling orifice and a supersonic free jet 
is formed ||T9ll . Several different effects are responsible for the change in composition 
||20|: radical recombination at external probe surface, acceleration into probe orifice, 
chemical relaxation in free jet, radial diffusion in free jet, skimmer interference and 
mach-number focusing. Probably, one of the most important effects is radial diffusion 
in the supersonic free jet. Huge radial pressure gradients exist just behind the 
sampling orifice which results into larger diffusion fluxes in the radial direction. 
The extent of this diffusion is different for different species, depending mainly on 
their mass, collision cross section and main collisional partner. Hence, pressure 
diffusion leads to mass separation downstream from the orifice. It has been shown 
that heavy species tend to remain on the central streamline because their diffusion is 
slow, whereas light species diffuse faster outwards along the pressure gradients. This 
effect leads to discrimination of light species against heavy ones in the composition 
of the MB. 

The influence of composition distortion is difficult to determine theoretically because 
it consists of many different effects. However, we can demonstrate it experimentally 
and use the observed changes in the signal to compensate for the composition 
distortion. 

The composition distortion is demonstrated by measuring the MS signals of a small 
amount of Ne (1 %) or N2O (0.5 %) in a He/air gas mixture. The air concentration in 
He/air mixture has been varied between pure He (light atoms with small collision 
cross section as main collision partners) and pure air (heavy diatomic molecules with 
larger collision cross section), whereas the Ne and N2O concentrations were kept 
constant. Ne is in pure He the heavier species and will enrich on the center line of the 
MB. While changing the He/air ratio towards pure air, Ne becomes the lighter species 
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Figure 4. Measurement of the Ne and N2O signal of different He/air mixtures with 
constant addition of Ne (1%) and N2O (0.5%) 



compared to the mean molecular mass of the mixture and will consequently diffuse 
faster out of the MB. Figure |4] shows the relative Ne signal measured by the MBMS as 
function of the air content in the He/air mixture. The Ne signal decreases by a factor 
of 4 during the He/air ratio variation. The N2O signal exhibits surprisingly a very 
similar behavior, even if it is more than twice heavier than Ne and heavier as N2 and 
O2 molecules in air. The difference between the masses of Ne, N2O and air (N2, 02) is 
probably not significant and the effect is dominated by their difference to He. 
The measurements in the figure S] show clearly that the composition distortion has to 
be taken into account. The amount of air in the effluent and the effect of composition 
distortion will strongly depend on the distance from the jet nozzle. The following 
procedure is adopted to take this effect into account in the signal calibration. First, 
the amount of air in the effluent at a given distance from the jet is determined by 
measuring the He and N2/O2 signals. These signals are compared to the measurements 
obtained with He/air mixtures with well defined composition. The local air content 
in the effluent can be determined in this way. This air content, together with the data 
in the figure IH is used to determine the correction factors necessary to compensate 
for the difference in the composition distortion in pure He and in He/air mixtures. 
The O and O3 signals are hence scaled up to the values, which would be obtained 
in pure He. The additional calibration is performed as described previously for the 
measurements in pure He |[T5|. Just briefly, O and O3 is calibrated by comparing their 
signals measured at masses 16 and 44 to signals measured for known densities of Ne 
and N2O in He, respectively. These gases have been chosen because of their similar 
mass and structure in comparison to O and O3. The electron impact ionization cross 
sections and the known mass dependence of the transmission function of the MS for 
different species are taken into account as well. Atomic oxygen has been measured 
with an electron energy in the ionizer of about 15 eV to prevent dissociative ionization 
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of molecular oxygen. Ozone has been measured with an electron energy of 70 eV 
because the signal does not overlap with dissociative ionization of other species. 



3. Fluid model of plasma effluent 

Absolute atomic oxygen and ozone densities in the effluent of the /i-APPJ obtained by 
MBMS can be used to reveal possible reaction pathways leading to their formation. 
The densities of ROS in the effluent of similar plasmas have already been measured 
and the chemical kinetics of these species have been modeled in the past | |2Tll22l . 
The reaction pathways proposed in these works are used here as a starting reaction 
scheme used in a 2D axially symmetric fluid model of our /j-APPJ. The model is 
based on the model already described in our previous work ||T3l , however, the effect 
of ambient air is added now and the measured absolute densities of O and O3 are used 
to validate the model results. The 2D axially symmetric geometry is used to be able 
to model the diffusion of air into the gas stream emanating from the plasma jet and 
to take the diffusion of ROS into account. The fluid model combines the solution of 
the Navier-Stokes equations for the gas flow in the gas mixture of He and air with the 
model of chemical kinetics of ROS. The model is solved using commercial COMSOL 
3.5 software. 



3.1. Model of the gas flow 

The gas flow through the jet and in the effluent is solved first. Only He and air 
(represented by N2, see later) are considered in the model of the gas flow, because the 
concentration of O2, which is injected into the plasma, and all other plasma products 
is less than 1 % and can be neglected. The gas flow is described by incompressible 
momentum conservation and continuity equations: 



- V • r](Vu + (Vu)^) + pu + Vp = (1) 
V • u = (2) 

with p being gas density, u gas flow velocity, i] dynamic viscosity and p the pressure 
(101325 Pa). No volume force is assumed to work on the gas. 

The transport of He through the ambient air is simulated as a diffusion-convection 
transport: 

j^+V- (-Dh.Vc) = -uVc (3) 

with c being the He concentration in the gas mixture (c G [0, 1]) and Dne being the 
diffusion coefficient of He in the gas mixture. The gas density, dynamic viscosity and 
also the diffusion coefficient are a function of the He concentration, which varies in 
space due to air diffusion from the sides into the effluent. The viscosity, density and 
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diffusion coefficients of N2 instead of air are used in the flow simulation for simplicity. 
The density of the He-N2 gas mixture is calculated as: 

pmixtiire = pHe ' C + pN^ ' 0- - c) (4) 

where c is the He concentration, p^e = 0.164 kg/m^ and Pn2 = 1.146 kg/m'^. The viscosity 
of the He-N2 mixture measured at 303.15 K is taken from the literature ||23l and fitted 
with following function: 

7](c) = -7.912xl0-^-c^+1.154xl0-^-c3-4.906xl0-'^-c2+3.355xl0-^-c+1.8xl0- 

The diffusion coefficients of He in the He-N2 gas mixture is calculated using following 
formula: 



Dfie in mixture ^He in He ^He in N2 

A self diffusion coefficient of He in He DneinHe = 1.7513x10"^ m^s"^ and a diffusion 
coefficient of He in N2 DneinNi = 0.7337x10"^ m^s'^ is used H241I251I. 
The jet and its effluent is modeled as two cylinders representing plasma channel and 
the effluent region as indicated in figure |5l The upper cylinder corresponds to the 
last 10 mm of the plasma channel and has a radius r=0.564 mm (1 mm^ cross-section 
area). The top base serves as a gas inlet, where a parabolic velocity profile through 
cylindrical tube calculated from the gas flow of 1.4 slm through 1 mm^ area (average 
velocity ~25.6 m/s) and He concentration equal one are used as boundary conditions. 
The bottom base serves as a gas outlet from the plasma channel into the effluent. No 
slip for the gas flow and insulation/symmetry for the He transport are selected as a 
boundary conditions for the outer wall of this cylinder. 

The effluent is modeled as a second cylinder with 3 mm radius and length of 50 mm. 
The bottom base and the outer part of the top base of this cylinder are assumed to 
be solid with no slip boundary condition. The outer wall is divided into two parts. 
The lower half of its length is defined as a gas outlet (boundary condition: convective 
flux, constant pressure p = 101325 Pa) and the upper part is defined as a gas inlet with 
constant gas velocity of 0.2 m/s (corresponding to a constant flow of 56 seem). He (c=l) 
or N2 (c=0) are introduced into the effluent region through this boundary and allow 
to simulate different ambient atmospheres. The inflow of ambient gas represents 
the "infinitely" larger reservoir of ambient air in our laboratory compared to the 
small volume of the plasma jet. Additionally, the flow field in the effluent region is 
stabilized by this inflow without formation of any vortexes in the numerical solution, 
which is a necessary condition to obtain a converging solution. This additional gas 
inflow into the effluent region could potentially influence the model results, but we 
have checked that selecting other gas velocities (4 time lower or higher) do not affect 
the model results presented later. 

The gas flow streamlines and a color map of the He concentration c for the case of 
air (N2) as ambient gas are also shown in figure |5l It can bee seen that air diffuses 
into the helium effluent as the distance from the nozzle of the /j-APPJ is increased. 
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Figure 5. Geometry used in the model with He concentration c in air and gas flow 
streamlines (white). Note the different scaling of r- and z-axis. 

Additionally, the results show that back diffusion of air into the plasma channel can 
be neglected in this modeled ideal case. 



3.2. Model of chemical kinetics 

When the gas flow and He concentration simulation is finished, the reactions and 
transport of ROS can be simulated. 

Again, the transport of given species z {i = 0, O3, 02(^Ag), or vibrationally excited 
oxygen 02(u)) is simulated as a diffusion-convection transport: 

-^ + V-{-D,yn,) = Ri-uVni (7) 

with Di and Ri being diffusion coefficient and production/loss term of a given 
particle due to gas phase reactions, respectively. Diffusion coefficients of O and 
O3 in He are Doi„He = 1.29x10"^ m^s"^ and Do,i„He = 0.713x10"^ m^s'^ IH. Do,i„N2^ 
0.137x10"'^ m^s~^ is the diffusion coefficient for O3 in N2 and has been averaged over 
the two stated values by Massman ||27||. Since no diffusion coefficient for O in N2 
could be found in literature, the diffusion coefficient for Ne in N2 has been calculated 
using the results of Chapman and Enskog 1125 1 . The masses of Ne and O are similar 
and both are atoms. This supports the assumption that both atoms exhibit a similar 
diffusion coefficient. A diffusion coefficient for Ne in N2 of DwemNz =0.338x10"'* m^s"* 
has been derived and used for the diffusion of O in N2. A gas temperature of 300 K 
is assumed in the model. The values of the diffusion coefficients of O2 in He {Do2inHe 
= 5.51x10"^ m^s"^) and in N2 (Do^mNz = 2.32x10"^ m^s"*) are used for 02(*A^) and 
vibrationally excited oxygen Oiiv). The dependence of D, on the He concentration 
c is calculated in the same way as described in the previous section (formula [H]). 
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reactions 



rate constants 



ref 



Ri O + O2 + He 

R2 O + O2 + O2 

R3 O + O2 + N2 

R4 O + O + He 

R5 + + O2 

Re O + O + N2 

R7 O + O3 ^ 

Rs 02CAg) + 03 



^ O3 + He 
^03 + 02 
^ O3 + N2 

O2 + He 
-^02 + 02 

O2 + N2 
O2 + O2 
-^0 + 202 



3.4 X 10-^'^(300/T^)i-2 m^s-i 
6 X 10-4^(300/r^)2-« m^s-i 
6 X 10-4^(300/T^)2-^ m^s-i 
1 X 10"^^ m^-i 
2.56 X 10-^^{300/Tgf^^ m^s"^ 
2.56 X 10-4'^(300/r^)0-^3 j^^g-i 

1.5 X 10-17 exp(-2250/r^)m3s-i 
6.01 X 10-17 exp(-2853/rg)m3s-i 



El 
est. ® 



in 

est. ® 



[3211 



Table 1. Reactions of Model 1 with rate constants. (® No literature values for three body 
reactions with N2 have been found. Instead, the reaction rate of the corresponding 
reaction with O2 as collision partner has been chosen, cf. ||3T1 .) 





reactions 


rate or time constants 


R9 


02{v) + O2 ^ O3 + 


2 X 10-21 m^s-i 


Rio 


02{V) ^ O2 


5 X 10-^ s 



Table 2. Additional reactions of Model 2. Rate and time constants have been used as 
fitting parameters. 



The dissociation degree of O2 molecules in the plasma is less than 5% under our 
experimental conditions. Therefore, the O2 density is not calculated in the model, but 
it is assumed to be constant (concentration of 0.6%) in the case of He as ambient gas 
and is calculated as noj = [0.006 + (1 - c) x 0.2] x no with no = 2.45 x 10^^ m-^ in the 
case of air as ambient gas. 

The boundary conditions are selected as follows. The boundary condition at any 
wall is different for stable O3 (insulation/symmetry) and reactive or excited species 
(O, 02(1 A^), 02(u)). Oxygen atoms and 02(u) can recombine or be deexcited at the 
surface, which has to be considered in the model. The Neumann's boundary condition 
considering the surface loss is therefore used lfT3ll28| . A surface loss probability of jS 
= lO-"' has been chosen for atomic oxygen and for 02{v). Due to the long lifetime of 
□2(1 Ag) ||29B , it has been assumed to be unreactive at the walls. Particles are introduced 
into the model by setting their density at the gas inlet into the plasma channel. They 
leave the volume through the convection losses at the gas outlet in the effluent region. 
Tabled] contains the gas phase reactions, which are implemented in the model. These 
reactions have been proposed by Jeong et al.flll and we use them in the simulation 
called Model 1. No reactions with ions or electrons are assumed because their 
densities in the effluent are negligible pO|. Additionally, the reactions shown in table 
|2]have been added to explain our measurements with air as ambient atmosphere. The 
reactions R9 - Rio of Model 2 are motivated and explained later in the chapter Results. 
The model with the complete set of reactions from tabled] and |2] is called Model 2. 
The volume is again divided into two parts: the plasma channel, where steady 
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value 



note 



O2 content in He/air mixture 



initial OiC^A^) density 



initial O density 
initial O3 density 
initial 02{v) density 



initial He density 
initial O2 density 



2.48 X 10^5 m-3 

1.49x10^3 
0.2 • (1-c) 

8.5xl020m-3 

3.5xl020m-3 
3 X 10^1 m-3 
7xl020m-3 



fitting parameter 
evaluated using |jT5l] and II32II 



extrapolated value from | [T5| 



extrapolated value from IITSl 



He at 101325 Pa 
0.6 % of He density 
20 % O2 in air 



gas temperature T, 



300 K 



Table 3. Used parameters in the model 



state situation with constant densities of all species is assumed, and effluent, where 
the reactive species can recombine. This assumption is supported by two-photon 
absorption laser-induced fluorescence measurements |[T4||, which show a constant O 
atom density in the plasma channel and a decaying density in the effluent. Therefore, 
no reactions are assumed in the plasma channel (the ROS densities are introduced into 
the model via boundary condition). The reactions from table [1] and |2] are considered 
only in they effluent region. 

Table |3] shows the used boundary conditions and parameters of both models. The 
shown initial densities are introduced at the gas inlet of the plasma channel. 



The experimentally determined and simulated O and O3 densities will be presented, 
compared and discussed here. First, the measurements with He as ambient gas will 
be shown and discussed. Afterwards, the results with ambient air will be presented. 
The results of Model 1 are discussed together with these measurements. Finally, the 
experiments performed to reveal the O3 generation in the effluent will be discussed. 
Here, the reactions of Model 2 are also motivated and explained and the corresponding 
calculated densities are discussed. 

4.1. Validation of the model with He as ambient atmosphere 

Figure [6] shows the atomic oxygen and ozone density measured by MBMS and two- 
photon absorption laser induced fluorescence spectroscopy (TALIF) in a helium 
ambiance as function of the distance. These measurements have already been 
discussed within a previous publication |T5| . In summary, the maximum O density 
is in the range of 7x10^^ cm"'^ - 9x10^^ cm"'^ and declines from its maximum with the 
distance. The O3 density is 3.7x10^^ cm"^ at the nozzle of the /j-APPJ and increases 
continuously with increasing distance. A saturation value of 1.4x10^^ cm"^ is reached 
beyond a distance of 40 mm. The atomic oxygen results obtained by MBMS and TALIF 



4. Results 
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5 10 15 20 25 30 35 40 45 50 

distance [mm] 

Figure 6. Atomic oxygen (measured by MBMS: •, measured by TALIF: O) and ozone 
(■) density as a function of the distance from the jet in ambient helium (applied 
electrode voltage: 230'Vrms, gas flow: 1.4 slm He with 0.6% O2). The lines reflect the 
densities derived by the Model 1 and 2. The TALIF measurement has been scaled 
down by a factor of 0.27. The results from MBMS and TALIF measurements are 
already published flSl . 

are in good agreement with each other despite of a discrepancy by a factor of 3.7. 
The TALIF measurements reveal a O densities 3.7 times higher than the one obtained 
by MBMS. The origin of the discrepancy are possibly uncertainties of the respective 
calibration processes. However, considering the fact that both measurements are 
performed by two different, independent diagnostics, this result can still be regarded 
as a very good agreement. 

The atomic oxygen density simulated by the fluid model (Model 1) matches perfectly 
the MBMS and TALIF measurements. In contrast, the simulated ozone density profile 
("Model 1 without filling effect" in figure[6]) doesn't correspond to the measurements. 
The modeled density exhibits a maximum in a distance of 21 mm and decreases 
beyond. For this model, helium is injected through the gas inlet at the side. This 
condition reflects a situation with the plasma effluent located in an infinitely large 
volume filled with He. However, the /j-APPJ is located in a chamber for the controlled 
He atmosphere with a volume of 1.1 dm^. This chamber fills continuously with 
ozone during the time the plasma is operated and the ozone diffuses back into the 
plasma effluent and causes an artificially higher ozone signal. To take this filling 
of the chamber with ozone into account, we assume that all O recombtnes into O3 
(neglecting O losses due to reactions R4-R8) and that an O3 density of no,i„it + t^OsMU 
is recycled through the gas inlet in the effluent region (boundary condition at the 
gas inlet). This assumption is reasonable, because we have tested the effect of the 
reactions R4-R8 on the model results under our experimental conditions. The change 
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distance [mm] 

Figure 7. Composition of the effluent as a function of the distance (1.4 slm helium 
flow) 

of the densities is less than 5% when reaction R4-R8 are implemented. The gas phase 
chemistry is dominated by the O recombination with O2 forming O3 (R1-R3). The 
result of the simulation with this recirculation of ozone is represented by the solid 
line in figured Now, the modeled data reproduce the trend of the measurement very 
well, even if the absolute values do not match perfectly. The profile of the O density 
is the same for both cases. 

In summary, we state that Model 1 and the measurements are in reasonable agreement 
and that the proposed chemistry of Jeong et al. can indeed reasonably well describe 
the O and O3 densities in the effluent in a helium ambiance. Here, the reactions Ri - R3 
would even be sufficient to simulate the effluent in an ambient helium atmosphere. 

4.2. Validation of the model with air as ambient atmosphere 

A measurement of N2 and O2 signals as function of the distance has been performed 
in ambient air to determine how fast air mixes into the helium effluent. The /.i-APPJ 
was just operated at a gas flow of 1.4 slm helium without addition of oxygen and 
without igniting a plasma. The distance between the ft-APPJ and the sampling orifice 
was varied between 1 and 52 mm while recording the He-, N2- and O2 signal (cf. 
figure [Z|. The sum of the N2 and O2 signal determines the air signal. The shown 
results have been corrected for composition distortion by measurements of different 
known He/air mixtures. 

Almost no air admixes into the helium effluent up to a distance of 4 mm. Air appears in 
the helium effluent beyond this point and its concentration is increasing continuously 
as the distance is increased. The gas composition on the symmetry axis, as calculated 
in the fluid model, is also shown in figure [71 The results from the model are in a very 
good agreement with the measured data and corroborates that the fluid model can 
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distance [mm] 

Figure 8. Atomic oxygen (•) and ozone (■) density as a function of the distance from 
the jet in ambient air (appHed electrode voltage: 230 V rms/ gas flow: 1.4 slm He with 
0.6 % O2). The solid lines reflect the densities derived by the model. 

very well simulate the admixture of air into the plasma effluent. 
4.3. Effluent chemistry in ambient air 

The MBMS measurements of the plasma effluent are now repeated with air as ambient 
gas. Again, the densities of atomic oxygen and ozone are determined and are shown 
in figurelH The atomic oxygen density is about 8x10^^ cm"'^ at the nozzle of the /.i-APPJ 
and declines with increasing distance almost linearly. But still in a distance of 15 mm 
oxygen atoms are present in the effluent (ca. 1.5x10^^ cm"'^). The atomic oxygen 
trends obtained in ambient air and helium are very similar. No TALIF measurements 
of O density are available for comparison because the TALIF signal is influenced by 
additional non radiative quenching of the excited states by air. 

There is no difference between the ozone densities as measured in He and air 
atmosphere up to a distance of 4 mm. This is an expected result because the air 
starts to admix into the effluent only beyond this distance (cf. figure At larger 
distances, the ozone density in ambient air significantly increases in contrast to the 
density in ambient helium. A maximum O3 density of 3.4xl0^^cm"'' is reached in a 
distance of 20 mm. Afterwards, the O3 density decreases to a value of 0.7xl0^^cm"^ 
at a distance of 48 mm. The steep decline of the O3 density in ambient air beyond 
a distance of 20 mm can be explained by the effect of buoyancy force. Absorption 
spectroscopy measurements of ozone by H. Bahre et al. ||33l have shown that the 
helium effluent (including the small concentration of ozone molecules) starts to 
ascend in a surrounding air atmosphere beyond a distance of ca. 20 mm. The origin 
of buoyancy can be the lighter density of helium in contrast to air combined with 
the slightly higher temperature of the plasma effluent. Later measurements of ozone 
with the modified MBMS setup does not show such a distinct shape due to a different 
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configuration. The jet is mounted vertically and the front plate with the sampling 
orifice is oriented horizontally at the modified MBMS setup. 

We have seen in the previous case, that the conversion of all O into O3 leads to a ozone 
density of about 1.2 to 1.4x10^^ cm"''. This means that more ozone is produced due 
to admixture of air into the effluent than due to atomic oxygen emanating from the 
plasma. However, there is no reaction in the proposed scheme of Jeong et al, which 
could explain the additional production of ozone. This can be clearly seen, when the 
measured densities are compared to the results of Model 1, also plotted in figure [8l 
The Model 1 predicts a faster decrease of the atomic oxygen density due to a faster 
reaction rate of reactions Ri, R2, R3, R5, and Re. This is caused by more available O2 
and N2 in the effluent because of the air diffusion into it. Correspondingly, ozone is 
produced faster compare to the case in ambient He. However, its maximal density 
predicted by Model 1 is still limited by the limited amount of oxygen atoms from the 
plasma and cannot exceed the limit given by the sum of initial O and O3 densities. 
Both O and O3 densities simulated by Model 1 do not reproduced the measured data 
(cf. figure ID). The atomic oxygen density does not decrease as fast as expected by the 
model, but stays relatively high even at a larger distance from the jet and, as already 
mentioned, much more ozone is produced. These results indicate that some reactions, 
which generate additional O atoms and O3 molecules, are missing in the simulation 
Model 1. 

The additional production of O atoms has already been observed by other authors. 
Reuter et al. | l22j l have observed that O atoms could be detected by TALIF even at the 
distance of 10 cm from a similar, but larger, planar jet operated in He with 0.5 % O2 
gas mixture and in He as ambient atmosphere. They concluded that electrons and 
ions as well as metastable helium atoms can be excluded as an origin for an energy 
transport into the effluent because their respective densities quickly decrease, as 
soon as electric excitation ceases. Additionally, they have detected energetic vacuum 
ultraviolet (VUV) radiation by OES and have shown that (V)UV radiation originating 
from the discharge region reaches far into the effluent. Therefore, they have concluded 
that "(V)UV radiation produces atomic oxygen in the effluent by dissociation of ozone 
or molecular oxygen." |2ZI This conclusion has been corroborated in their following 
work, where the plasma emitted photons have been blocked, which resulted in a 
decrease of O densities in the effluent Il34l . (V)UV radiation as a source of reactive 
species in the effluent of the jU-APPJ has also been observed indirectly during the 
inactivation of bacteria by the effluent of the X-Jet ||T3l[T8l . We will try now, with the 
help of the absolute calibrated data, additional MBMS measurements, and adjustment 
of the fluid model, to reveal more details about a possible O and O3 generation scheme 
in the plasma effluent. 
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4.4. Revealing the origin of the additional O and ozone 
The possible effect of air back diffusion 

As already discussed, the results in the literature indicate that plasma generated 
(V)UV photons can initiate photodissociation reactions leading to formation of 
additional O and subsequently O3. Before we, however, discuss the possible 
photodissociation reactions, we need to check for the possible effect of air back 
diffusion into the plasma channel. The gas flow simulation has shown, that air 
back diffusion is almost negligible. But the model represents an ideal case without 
effects due to actual design and realization of the plasma jet such as diffusion along 
the corner edge of the square geometry. Additionally, the simulation shows that 
even in the ideal case the air can come into direct contact with the active plasma 
zone. Point A in figured marks the edge of the plasma channel, where it takes place. 
The electron driven dissociation of atmospheric O2 followed by the diffusion of the 
products towards the jet axis could also be a source of additional O and O3. To exclude 
this effect, we have built a modified jU-APPJ with an elongated gas channel. This gas 
channel is made from glass plates, starts directly at the electrodes, has a length of 
3 mm and a 1x1 mm^ cross section. The plasma effluent is therefore confined the 
first 3 mm after the plasma region and a direct interaction between atmospheric air 
and the plasma is prevented. Figure |9] shows the comparison of the ozone signals 
obtained with the standard /j-APPJ without extension and the extended /i-APPJ. Only 
relative data are shown, no absolute calibration and compensation for the changing 
composition distortion have been performed. No significant difference can be found 
between both measurements. One can, therefore, conclude that diffusion of air into 
the plasma cannot explain the increased ozone (and also O) production. 
During these measurements, the newer MBMS setup with horizontally oriented front 
plate has been used. Therefore, the influence of the helium buoyancy is different 
and the shape of the measured ozone signal differs a little bit from the previous 
measurement (cf. figure H]). 

Photodissociation as the source of additional O and O3 
Photodissociation of oxygen molecules or ozone has been suggested as a possible 
source of additional O and consequently also O3. The quantitative and space resolved 
MBMS measurements of O3 allows now the quantitative estimation of the necessary 
O3 production rate needed to explain the rise of its density as observed in figure [8l 
This estimated production rate can be used further to estimated the necessary photon 
flux, if the photodissociation is the main source of this ozone. The O3 density rises 
from 0.5x10^^ to 3.0x10^^ cm"'^ on the distance of 13 mm (starting at 6 mm from the jet 
nozzle). The fluid model predicts the average gas velocity on the axis of about 25 m/s 
resulting in the transport time t=0.013/25 s = 0.52 ms. The averaged ozone production 
rate in the effluent can therefore be estimated as: 

= 2.5x10^^ cm-3/5.2xl0-^ s = 4.8x10^^ cm" ~ 2-T ■ a ph-n source -(8) 
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Figure 9. Relative MS signal at mass 48 amu (ozone) as function of the distance to 
the nozzle for the standard |ti-APPJ and extended /J-APPJ under same conditions (in 
ambient air, applied electrode voltage: 230 VrmS/ gas flow: 1.4 slm He with 0.6% O2). 



The right hand side of equation|8]is the estimated production rate of the ozone due to 
photodissociation of some source molecule (e.g. O2 or O3) with density nsource- The Oph 
is the photodissociation cross section, T is the photon flux from plasma per unit area 
and second, and factor 2 takes into account that two O3 molecules can be produced 
in each case. The possible reactions could be the direct dissociation of O2 into 2 O 
atoms, which than quickly react with O2 from the air into O3 or the photodissociation 
of O3 into O and highly vibrationally excited 02{v), which can again react with O2 
producing O3 and O in this case. The two O atoms will as in the previous case react 
to O3 providing again a net production of two new O3 molecules. 
The typical maximum of the photodissociation cross sections of O2 or O3 molecules 
(Oph ~ 10"^^ cm^) p5l l36ll , can be used to estimate the lowest necessary photon flux 
in the ideal case, in which each photon is absorbed with maximum probability. With 
no2 ~ lO^'^m-^andnog ~ 2x10^^ m"Hhe photon fluxes of at least Tq^ ~ 2.5x10^^ m'^s"^ 
and Toj ~ 1.2 X 10^^m"^s"^ are necessary for the ozone production via photodisso- 
ciation of O2 or O3, respectively. With the 1 mm^ cross-section area of the jet and a 
chosen photon energy of 5 eV, we obtain a necessary energy flux of at least 2 mW in 
the case of dissociation of O2 and 1 W in the case of photodissociation of O3. The 
latter reaction can be directly excluded, because the power absorbed by the plasma 
is for sure lower than 1 W and the plasma cannot generate so many photons. Even 
in the case of photodissociation of O2 seems the power of the photon flux too high, 
considering the crude simplifications in the estimation of the photon flux. However, 
it could still contribute to the production of O and ozone in the effluent. 

Luckily, the photodissociation reaction of O2 by plasma generated photons can di- 
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rectly be tested by the X-Jet modification of the fi-APPJ. For this experiment the 
additional He flow with different concentrations of O2 is flown through the side 
channel of the X-Jet, whereas the plasma is ignited in the standard mixture of He with 
0.6% of O2 in the direct channel between electrodes. The additional flow diverts the 
plasma effluent from the direct channel into the side channel after the crossing of both 
channels and fills at the same time the part of the direct channel after the crossing 
with a He/02 gas mixture as shown in the figure |2^. The plasma generated photons, 
which are propagating through the direct channel of the X-Jet, overlap spatially with 
the He/02 gas- The possible flow of O and O3 from the direct channel, generated now 
by photodissociation reactions, can be measured by the MBMS system. O2 concentra- 
tions of up to 20% and different distances between the nozzle of the direct channel and 
the sampling orifice of the MBMS were tested to allow the photons traveling different 
lengths through the He/02 mixture and also to allow air to diffuse into the gas flow. 
Ozone has never been detected in these measurements. Additionally, the normal 
/.i-APPJ operated in pure He gas, without addition of O2 gas, has been operated in 
ambient atmosphere and the production of ozone in the effluent at different distances 
from the jet has been measured. Even more VUV radiation is expected to be produced 
in this case due to slower quenching of excited states and higher densities of helium 
excimers. The ozone density was also in this case below the detection limit of the 
MS. All these measurements are evidence that photodissociation of ground state O2 
molecules can be very probably excluded as the source of additional O and O3 in the 
effluent. It seems, that photodissociation reactions are not the source of additional O 
atoms and ozone molecules under our experimental conditions and some other oxy- 
gen containing species generated in the plasma have to be involved in this process. 

Summarizing the facts discussed up to now, the following conditions have to be 
fulfilled in the involved reaction mechanism: (i) the VUV radiation can be probably 
excluded as a source of additional O and O3 due to the above discussed reasons; 
(ii) molecular ground state oxygen from air is required because the additional ozone 
signal only appears when air diffuses into the effluent; (iii) atomic oxygen is produced 
in this reaction because the O densities are at larger distances from the jet higher as 
predicted by Model 1 (cf. figur^; (iv) some reactive oxygen species from the plasma 
have to be involved in the reaction because the interaction of the effluent of pure He 
plasma with ambient air does not lead to formation of O and O3; and (v) the reaction 
has to be somehow limited by the O2 density because the effect is weak with 0.6% of 
O2 in the mixture (the measurements with He as ambient atmosphere) but becomes 
significant at higher O2 densities when air admixes into the effluent. 

The reaction of the highly vibrationally excited 02{v) molecules, which are gener- 
ated in the plasma and have a limited life time, with the ground state O2 is proposed 
here as a possible reaction mechanism. The vibrationally excited oxygen molecules 
can directly react with ground state oxygen molecules (reaction Rg) producing an 
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oxygen atom and an ozone molecule. Two O3 are therefore effectively formed per 
reaction. The reactions of highly excited 02(t^>26) molecules have been for example 
discussed in the literature and were proposed to solve the ozone deficit problem in the 
atmosphere ll37ll38B . Additionally, the 02(u) relax decay to its ground state (reaction 
Rio) is proposed to explain the small effect observed in the experiments in ambient 
He. This reaction scheme has been included to the simulation (Model 2). Since some 
of the needed parameters are unknown, they have been adjusted until Model 2 fits 
reasonably well the measured densities. The fitting parameters are reaction rate of 
reaction R9 and the lifetime and the initial density of the vibrationally excited oxygen 
molecules. The used values of the fitting parameters are shown in the tables |2] and |3l 
The results of the Model 2 are shown in figure [6] for the case of ambient helium and 
in figure [8] for the case of ambient air. Model 2 is in very good agreement with the 
measurements for both cases, ambient helium and ambient air. 

This result corroborates that a reasonable reaction mechanism without involvement 
of photons can explain under our experimental conditions the observed results. We 
stress here, that it is just a hypothesis. A validation by possible measurements of 
02iv) and by comparison of fitted rates and relaxation times with measured rates 
and relaxation times, if they exist, should be performed before accepting this reac- 
tion mechanism. Additionally, it should be checked whether this proposed reaction 
mechanism can explain the observations, which have lead to the conclusion that VUV 
and UV photons are involved in the effluent chemistry. 

5. Conclusion 

Molecular beam mass spectrometry is used for measurements of absolute densities 
of O and O3 in the effluent of a He/02 micro-scaled atmospheric pressure plasma jet 
operated in ambient air and ambient helium. The effect of the composition distortion 
in the molecular beam is studied and carefully considered in the calibration process. 
Additionally, the admixture of air into the plasma effluent has been quantified. A 
simple fluid model of the gas flow, admixture of the air into the effluent, and reaction 
kinetics of several reactive oxygen species is constructed and used to test different 
reaction mechanisms. 

It was observed that the admixture of air into the effluent leads to an additional 
production of oxygen atoms and a fast increase of ozone density. Surprisingly, two 
to three times more ozone is produced in the additional reactions of the plasma 
effluent with the ambient air than what is observed in ambient helium without big 
differences in the distance variation of atomic oxygen density. This observation can 
have big implication for the use of similar jets in the plasma medicine and plasma 
sterilization applications. The possible role of plasma generated photons has been 
carefully tested. It is shown that the photodissociation of O2 and O3 is not responsible 
for the observed effect and that the photon flux is very probably too low to have 
any effect. A reaction scheme involving the reaction of plasma produced highly 
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vibrationally excited oxygen molecules with a limited lifetime with ground state O2 
molecules is proposed as a possible explanation of the observed O and O3 densities. 
A very good agreement between measured and simulated densities is achieved. 
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